Relaxation and Coarsening Dynamics in Superconducting Arrays 
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We investigate the nonequilibrium coarsening dynamics in two-dimensional overdamped super- 
conducting arrays under zero external current, where ohmic dissipation occurs on junctions between 
superconducting islands through uniform resistance. The nonequilibrium relaxation of the unfrus- 
trated array and also of the fully frustrated array, quenched to low temperature ordered states or 
quasi-ordered ones, is dominated by characteristic features of coarsening processes via decay of point 
and line defects, respectively. In the case of unfrustrated arrays, it is argued that due to fmiteness of 
the friction constant for a vortex (in the limit of large spatial extent of the vortex), the typical length 
scale grows as i s ~ t 1 / 2 accompanied by the number of point vortices decaying as N v ~ 1/t. This is 
in contrast with the case that dominant dissipation occurs between each island and the substrate, 
where the friction constant diverges logarithmically and the length scale exhibits diffusive growth 
with a logarithmic correction term. We perform extensive numerical simulations, to obtain results 
in reasonable agreement. In the case of fully frustrated arrays, the domain growth of Ising-like chiral 
order exhibits the low-temperature behavior t q ~ t l ' Zq , with the growth exponent l/z q apparently 
showing a strong temperature dependence in the low-temperature limit. 

PACS numbers: 74.50+r, 67.40.Fd 



I. INTRODUCTION 

In a statistical system quenched from the random dis- 
ordered state to a low-temperature state below the tran- 
sition temperature, the average length scale t of ordered 
domains typically grows in time as a power law I ~ t 1 '" , 
where the growth exponent X/z depends on the dimen- 
sion of the space and of the relevant order parameter, in 
addition to the conserved or noncenserved nature of the 
latter in the relaxation dynamics .a Depending crucially 
on the dimension of the order parameter, characteris- 
tic topological defects such as point vortices or domain 
walls are generated in the initial disordered state and 
the annihilation of these defects gives the main mech- 
anism of coarsening and phase ordering in the system. 
The observed self-similarity of these coarsening systems 
at different time instants, is embodied in the so-called 
dynamic scaling hypothesis of the equal-time spatial cor- 
relation function of the order parameter. 

One of the extensively studied cases is the system with 
the order parameter dimension n = 2 in two spatial 
dimensions d = 2, i.e., the two-dimensional (2D) XY 
model, rwheije the relevant topological defect is the point 
vortexE&QCJ A convenient physical realization of the 2D 
XY model can be found in the 2D regular arrays of su- 
perconducting Josephson junctions, where the phases of 
the superconducting order parameters correspond to the 
XY spins. In such a superconducting array, frustration 
may he induced simply by applying an external magnetic 
field.B It affects thermodynamic properties ia a crucial 
way, leading to a variety of phase transitions!] 



In view of the remarkable diversity in the equilibrium 
properties, also expected for the system are a variety of 
interesting dynamic behaviors, particularly in relaxation 
toward equilibrium. The quasi-long-range order present 
in the 2D pure XY modelja describing the array with- 
out frustration (i.e., in the absence of a magnetic field), 
is expected to exhibit coarsening through annihilation of 
vortex-antivortex pairs, giving rise to characteristic alge- 
braic relaxation. The vortex dynamics and the resulting 
scaling behavior in the 2D XY model has been an issue 
for some time, concerning late-time scaling and true na- 
ture of the asymptotic growth law. It is also of interest 
to examine the effects of the additional long-range or- 
der on the relaxation in the fully frustrated array (with 
half the flux quantum per plaquette) J3'B where the relax- 
ation dynamics is expected to be dominated by coarsen- 
ing processes accompanying decay of domain walls and 
also point defects; the latter corresponds to the corners 
on the domain walls. 

The main purpose of this paper is to probe the proper- 
ties of the nonequilibrium relaxation dynamics associated 
with coarsening processes in the unfrustrated and fully 
frustrated superconducting arrays under rapid quench- 
ing from the high-temperature disordered state to the 
low-temperature ordered or quasi-ordered states. We 
would like to understand what kind of specific features 
emerge in the relaxation dynamics of these systems gov- 
erned by the resistively-shunted-junction (RSJ) dynamic 
equations, especially in connection with the coarsening 
processes. The time-dependent behaviors of various one- 
time and two-time correlations involving energy and or- 
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der parameters such as vortex density and chirality are 
investigated numerically. 

In the case of an unfrustrated array, we argue that the 
friction constant of a slow moving vortex remains finite 
in the limit of large extent of the vortex;E2l this leads (via 
a simple force-velocity relation) to a power-law growth 
of length scales with growth exponent 1/2, together with 
the excess number of vortices decaying as N v ~ t~ a " 
with a v — 1. Extensive simulations indeed show that 
a v takes a value that is very close to unity in a wide 
range of temperatures. This is in contrast to the simula- 
tion results of the ordinary phase ordering process based 
on either Monte Carlo algorithms or simple Langevin dy- 
namics of the XY model (see section III for details), 
where a logarithmic correction factor is observed in the 
time dependence of the growing length scale. The latter 
case corresponds to the situation in Josephson-junction 
arrays where dissipation occurs between each island and 
the substrate. 

In relation to the coarsening dynamics of fully frus- 
trated arrays, Lee, et alH3 performed dynamic simu- 
lations on the fully frustrated XY models, based on 
the simple Langevin equations for XY phase angles. 
It was shown that the domain growth exponents are 
temperature-dependent and that there exist two regimes 
of domain-growth morphology. Later, coarsening was 
also studied in the dual representation, i.e., the lattice 
Coulomb gas of fractional vortices, through the use p£ 
Monte Carlo methods on the vortex degrees of freedom.La 
The latter work, even though based on different dynam- 
ics, confirmed most features of the previous results and 
gave some numerical evidence for unbinding of the so- 
called single-step kinks at a finite temperature. Those 
rich features of coarsening in the fully frustrated XY 
model can be partly attributed to the existence of the 
Ising-like chiral degrees of freedom in addition to the un- 
derlying phase degrees of freedom. Here, in this work, we 
perform dynamic simulations on the coarsening of the 
fully frustrated Josephson-junction arrays, governed by 
the RSJ equations, and examine whether temperature- 
dependent coarsening also emerges. It is found that, simi- 
larly to the previous works, the domain growth exponents 
are strongly temperature-dependent: The exponents are 
proportional to the temperature in the low-temperature 
limit and increase monotonically with the temperature. 



where we have employed the abbreviations cf>ij = 
(jji—cjjj—Aij and = r^— and the primed summation 
runs over the nearest neighbors of grain i. The position 
of grain i is represented by = (xi,yi) with the lattice 
constant set equal to unity while the gauge field Ay is 
given by the line integral of the vector potential A: 
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with the flux quantum <I>o = hc/2e. The frustration pa- 
rameter /, which measures the number of flux quanta per 
plaquette, is given by the directional sum of the gauge 
field Ay around a plaquette: 
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In Eq. (|]J) the energy and the time have been expressed 
in units of HI c /2e and h/2eRI c , respectively, with single- 
junction critical current I c and shunt resistance R. The 
thermal noise current Qj is assumed to be white, satisfy- 



ing 



((ij(t+T)(ki(t)) = 2k B T5(T){8 %k 5 ]l -5 ll 5 ]k ) (4) 



at temperature T. Henceforth we set the Boltzmann 
constant fc# = 1. The dynamics of the twist variables 
A = (A X) A y ) is governed by 



dA a 1 ^ . 



A Q )+Ca, 



(5) 



where denotes the summation over all nearest- 

neighboring pairs in the a-direction (a = x, y) and £ a 
satisfies 



2T 

(Ca(t+r)( b (t)) = j^S ab 6(r). 



(6) 



Note that in equilibrium the above set of equations of 
motion leads to a Gibbs distribution with the Hamilto- 
nian for the frustrated XY model 



H = — cos 



r«-A), 



(7) 



II. EQUATIONS OF MOTION 

We begin with the set of equations of motion for the 
phases {4>i} of the superconducting order parameters in 
an L x L square array. In the RSJpHiodel with the fluc- 
tuating twist boundary conditions ,E3 they read: 



dt 



0, (1) 



where the fluctuating twist boundary conditions have 
been incorporated. In numerical simulations, we have 
integrated directly the equations of motions in Eqs. (|l|) 
and (|) via the modified Eulcr method with time steps 
At = 0.05. Simulations have been performed on square 
arrays with linear size L = 400 (unfrustrated arrays) and 
L = 128 (fully frustrated arrays). 

To study the relaxation and coarsening in the super- 
conducting arrays, we let the systems evolve from random 
initial configurations at given temperatures and measure 
the following quantities: 
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1. The excess amount of topological defects such as 
vortices and domain walls 



N(t) = N(t) -N(oo), 



(8) 



where N(t) is the total amount (number or length) 
at time t and N(oo) refers to the amount in equi- 
librium. [In the case of vortices N denotes the sum 
of the total number of vortices (plus charges) and 
that of antivortices (minus charges).] It is under- 
stood that the appropriate ensemble average over 
random initial configurations is to be taken. 



2. Excess energy relaxation defined to be 
25(f) = 25(f) - E(oo) 



(9) 



where the energy (per site) 25(f) at time t is given 
by 



25(f) 



L 2 ^ 



K*) 



A(f) 



(10) 



with the summation over all nearest-neighboring 
pairs and (■ • ■ ) denoting the ensemble average over 
random initial configurations. 

3. The equal-time spatial correlation function 

C(r,*) = (O*(r,t)O(0,t)> (11) 

of the order parameter 0(r, f) = e l ^ r,t \ 

According to the dynamic scaling hypothesis, the corre- 
lation function assumes the 2D scaling form 



C{r,t)=r-«g{r/l) 



(12) 



with the appropriate scaling function g(x), where I = £(f) 
is the typical length scale of ordered domains (at time t). 
While the exponent rj vanishes (rj = 0) in the usual non- 
critical quenching to long-range ordered states, it takes 
nonzero values in the case of quenching to critical states, 
depending on the temperature T, i.e., r\ = rj(T). 



III. COARSENING AND VORTEX DYNAMICS 
IN UNFRUSTRATED ARRAYS 

It is well known that the unfrustrated system (Ay = 
0), described by the XY model in equilibrium, is criti- 
cal below the BerezinskiiJCosterlitz-Thouless transition 
temperature T c (w 0.89)o Previous works on relaxation 
of the XY model dealt mostly with the time-dependent 
Ginzburg-Landau equations (either the soft-spin version 
or the hard-spin version). The hard-spin version, where 
the magnitude of each spin is fixed, is described by the 
simple Langevin equations for phases: 



SH 



dt 



(13) 



where H = j) cos (0i — 0i) i s the Hamiltonian and 

the noise Q satisfies (Ci(*+ T )Ci(*)) — 2TSijS(r). Experi- 
mentally, this simple Langevin dynamics corresponds to 
the situation where dissipation through the resistance be- 
tween superconducting islands and the substrate domi- 
nates over those through the junction resistance. 

In this model with no frustration, it was argued that 
the friction constant 7 of a slow moving vortex depends 
logarithmically on the size (spatial extent) r of the vor- 
tex, i.e., 7 ~ ln(r/ro), where the cutoff length scale ro 
may be taken as the size of the vortex core.Ej Combining 
this result with the coulomb force law of F(r) ~ — 
we get the relation fln(r/ro) ~ — r _1 , which leads to 
the conclusion that the time dependence of the vortex 
number exhibits logarithmic correction to the power law: 
N v ~ f _1 hif. The typical length scale, e.g., the se 



ration between vortices, is then expected to followt2lt3 
£ s ~ (i/lni) 1 / 2 . Numerical confirmation of ikis predic- 
tion has turned out rather tricky and subtle.liS 

Here, we apply a similar argument to the coupled RSJ 
dynamics of the superconducting arrays. Interestingly, 
the mobility of a vortex whose motion is governed by 
the RSJ dynamics is found not to vanish in the large- 
separation limit (r — > 00) but to remain finite. This can 
be easily seen from the following argument based on the 
equations of motion in Eq. (|l|) for the phase variables. 
Equation (pi) can be rewritten in the following formLLZI 



SH 



dt 



(14) 



where —Mij is the discrete Laplacian matrix and the 
Hamiltonian H is given by Eq. In order to obtain 
the mobility of an isolated vortex, we compute the energy 
dissipation for an isolated vortex moving with a constant 
velocity v along, say, the x direction. The corresponding 
phase configuration takes the form (f> v (x—vt,y), where 
(j) v (x,y) represents the phase at site i = (x,y) in the 
presence of an isolated (static) vortex. At low temper- 
atures, where the noise term may be neglected, we use 
Eq. ( |14| ) and write the energy dissipation in the formliS 



dE 
~dt 



E 



SH dfa 
Scj)i dt 



E^ 



dt dt 



(15) 



Inserting the vortex solution <p v (x—vt,y) in the above 
equation and working in the continuum notation, we get 
the following expression for the energy dissipation 



dE 2 



(16) 



where the friction constant (or the inverse of the mobil- 
ity) j v is given by 



dxdyd x (j) v (x-vt, y)V 2 d x (f) v {x-vt, y). (17) 



In the derivation of Eq. ( |17| ) from Eq. (|15j), it has been 
noted that in the continuum limit, the discrete Laplacian 
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FIG. 1: Relaxation of (a) the vortex number N v at temper- 
ature T — 0.30 and 0.45 and of (b) the excess energy E at 
T = 0.45 in the unfrustrated array of size L = 400. In (a) the 
power-law behavior t -0 - 99 is also plotted whereas (b) displays 
the best fit to the form E — at^ 1 ln(_/_o) with a — 0.175 and 
t = 0.06. 



matrix —Mij turns into V 2 and the time derivative be- 
comes velocity multiplied by the spatial derivative with 
respect to x for the ansatz vortex configuration. In the 
limit of zero vortex velocity (v — > 0), one can easily see 
from dimensional argument that the mobility becomes fi- 
nite for a vortex of an infinite extent. Namely, the above 
integral for j v is convergent in the large-size limit, ow- 
ing to the presence of the additional Laplacian derivative 
term V 2 . Since we have d x cj) v (x,y) ~ r~ 4 for a single 
isolated vortex, the integrand in the above equation goes 
as ~ r~ 4 , which is infrared convergent in two dimensions. 

Since the friction constant converges quickly to a fi- 
nite value 70 as the extent of the vortex is increased, the 
equation of motion resulting from the force balance reads 



from which one can estimate the annihilation time for a 
vortex-antivortex pair of size £ to be t ~ (£/£o) 2 , with 
the cutoff length scale £o corresponding to the size of 
the vortex core. This can be interpreted as the growth 
law for the typical (average) separation between vortices 
as £ ~ t 1 / 2 , which in turn leads to the vortex number 
N v ~ £~ 2 n_ £ _1 . Figure |l](a) shows the relaxation of the 
number of vortices at temperatures T = 0.3 and 0.45. 
Observed in the late-time regime are good fits with N v ~ 
t~ " and t~ - 98 at T = 0.45 and at T = 0.3, respectively, 
in excellent agreement with the expected i _1 behavior. 

This result may be used to obtain the energy relaxation 
behavior: On the basis of the dynamic scaling hypothesis 
and the Porod lawj_ which refers to the power-law singu- 
larity in the short (scaled) distance limit of the equal-time 
spatial correlation function due to the presence of topo- 
logical defects, it is known that tie excess energy for the 
0{n) model with n = 2 scales asEJ 



E 



>•(} 



(19) 



where the typical length scale may be taken to be sep- 
aration between vortices and the cutoff to be 
the vortex core size (tq ~ £o) m °ur system. Noting 
the simulation result £ ~ t 1 / 2 presented above, we get 
E ~ t _1 ln(t/to) with io being approximately the ini- 
tial microscopic time scale taken by a vortex to move 
a single lattice spacing. The excess energy versus time 
measured at T = 0.45 is shown in Fig. |l|(b) and indeed 
fitted quite well to the expected form E as at^ 1 ln(_/_o) 
with a = 0.17 and io = 0.06. We have also tried power- 
law fitting of the form E ~ t~ a , to obtain the effective 
exponent a w 0.89. Such a value which is smaller than 
unity can be understood as coming from the logarithmic 
factor. 

Another test on the growing length scale in the unfrus- 
trated array can be performed by calculating the equal- 
time phase correlation function 



i+r 



(t)-0i(i)]>. (20) 



7o r 



Figure shows the equal-time spatial correlation func- 
tion of phases at temperature T — 0.45. Behavior of 
the phase correlation function C s (r, t) with distance r is 
displayed in Fig || at various time stages t = 8 to 2048. 
Critical dynamic scaling has been attempted with vari- 
ous values of r/: For each trial value of rj, the rescaled 
function C s (r,t) = r ?, C s (r, t) has been plotted in terms 
of the rescaled distance r/£ s , where the length scale t s 
as a function of time t has been obtained from the con- 
dition C s (r=£ s ,t) = Co = 0.4. The value of r/ yielding 
the best scaling collapse has been chosen, which is shown 
in Fig. ||(b). The resulting optimal value 77 = 0.083 is 
in fair agreement with the value from the theoretically, 
predicted form of the equilibrium correlation exponenttfj 
rj eq w T/2tt + T 2 /4tt w 0.0877. The length scale l s ob- 
(18) tained in this way is presented in Fig. 0(c), where one 
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FIG. 2: (a) Spatial behavior of the equal-time phase correla- 
tion function at various time stages in the unfrustrated array 
of size L = 400 at temperature T = 0.45. (b) Scaling collapse 
of the data in (a), plotted in terms of the rescaled correlation 
function with rj = 0.083 and the scaling length £ a defined in 
the text, (c) Scaling length versus time at T = 0.45, exhibit- 
ing power-law growth t s ~ t 0AS . 



finds that £ s ~ t l / z - with l/z s w 0.48, slightly below the 
expected value 1/2. We have also performed simulations 
at various temperatures, and found that the value of l/z s 
is nearly constant, showing only small variations in the 
range 1/ z s sw 0.47 to 0.48. The reason behind this small 
but apparently systematic deviation is not clear at this 
stage. Presumably, it is related to the fact that we are 
working with a lattice system where small local barriers 
exist for the motion of point vortices. At the early-time 
stage when the vortices are close together, these barriers 
are negligible compared with inter-vortex forces; at the 
late-time stage when vortices are far away, local energy 
barriers are likely to hinder appreciably the vortex mo- 
tion, possibly resulting in smaller values of the effective 
growth exponent. 



IV. COARSENING IN THE RELAXATION 
DYNAMICS OF FULLY FRUSTRATED ARRAYS 

We now consider the fully frustrated array (/=l/2) 
with the gauge field: 





TTXi 



for Tj 
for r, 



x 

•y, 



(21) 



which corresponds to the system of XY spins with one 
antifcrromagnetic coupling per plaquette. The fully frus- 
trated array has the discrete Z2 symmetry in addition 
to the underlying U(l) symmetry, yielding the antifcr- 
romagnetic long-range ocder in the chirality at temper- 
atures T < T c (« 0.45) .El Due to the existence of the 
chirality order parameter, the coarsening dynamics in a 
fully frustrated array proceeds in a manner quite different 
from that of the unfrustrated array. 

In order to probe nonequilibrium coarsening of a fully 
frustrated array, we measure the equal-time spatial corre- 
lations of the phase and the chirality. The phase correla- 
tion function in the fully frustrated array can be defined 
in the same way as Eq. (|2Q|) for the unfrustrated array, 
only if it is assumed that the correlations are measured 
only for even lattice spacings along both the axes. Re- 
stricting the position vectors to those with components of 
even lattice spacings corresponds to forming a superlat- 
tice, each plaquette of which consists of four plaquettes 
of the original lattice. Such a superlattice system yields 
a ferromagnetic ground state within itself at zero tem- 
perature. Meanwhile, the chirality correlation function 
is defined as follows: 



C q (r,t) 



R 



<?R+r(t)<7R.(i)), 



where the chirality 

qn(t) = sgn^sin 



/(*) 



A(i) 



(22) 



(23) 



with denoting the directional plaquette summation 
of links around dual lattice site R, describes the degrees 
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C q (r,t) 





FIG. 3: (a) Spatial behavior of the equal-time chirality correlation function at various time stages in the fully frustrated array 
of size L = 128 at temperature T = 0.15. (b) Scaling collapse of the data in (a) with the appropriate scaling length t q . (c) 
Scaling length versus time at various temperatures T = 0.3, 0.2, 0.15, 0.1, 0.05 from above, revealing temperature-dependent 
growth, (d) Growth exponent versus temperature, manifesting that growth becomes slower as the temperature is lowered. 



of freedom associated with the Z2 symmetry. The factor 
(— l) x+y , in the definition of C q (r,t), where r = (x,y), 
takes care of the staggered nature of the chirality order- 
ing. Due to the existence of these two kinds of degrees 
of freedom, the relaxation and ordering process proceeds 
via annihilation of two kinds of defects i.e., point defects 
(tending to restore the U(l) symmetry) and line defects 
(domain walls tending to restore the chiral Z2 symme- 
try). Point defects correspond to the corners of line dfe 
fects where fractional charges of magnitude 1/4 reside.ta 
Hence the annihilation processes of these two kinds of 
defects are closely coupled with each other. Note that 
conventional point vortices with integer charges, possess- 
ing high energies, are annihilated very quickly in the early 
stage of coarsening; thus they are irrelevant in the late- 
time coarsening behavior of the system. 

Here we present our simulation results for the ordering 
dynamics of only the chirality since the phase ordering 
exhibits similar features except for the additional r\ expo- 



nent due to the criticality. Figure g shows the equal-time 
correlation function C q (r, t) of the (staggered) chirality 
at temperature T = 0.15. Spatial behavior of C q (r,t) is 
exhibited in Fig. ||(a) at various time stages t = 16 to 
2048. Scaling of the correlations can be obtained simi- 
larly to the case of the unfrustrated array, except for the 
obvious fact that rj q = 0, reflecting the emergence of true 
long-range order. The resulting nice scaling collapse is 
shown in Fig. ||(b) . The scaling length l q is plotted ver- 
sus time in Fig. 0(c) at various temperatures, from which 
one obtains the relation l q ~ t 1 ^" 1 with the appropriate 
growth exponent l/z q . 

In contrast to the case of the unfrustrated array, the 
growth exponent in the fully frustrated array depends 
strongly on the temperature, as shown in Fig. [|(d). For 
example, we have the value l/z q w 0.38 at temperature 
T = 0.2 but it decreases to l/z q w 0.32 at T = 0.15. 
At low temperatures, in particular, l/z q depends almost 
linearly on the temperature, which has also been ob- 
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(b) (d) 

FIG. 4: Snapshots of the configuration of the staggered chiral- 
ity in the fully frustrated array of size L = 64 at temperature 
T and time t. (a) T = 0.1 and t = 32; (b) T = 0.1 and 
t = 1024; (c) T = 0.4 and f = 32; (d) T = 0.4 and t = 1024. 



served in the ordering dynamics of the fully frustrated 
XY model based on the simple Langevin equationO On 
the other hand, the temperature dependence of the do- 
main growth emerges only after some time (t >t q rj 30); 
before that time (t < t q ) the domain growth apparently 
does not change with the temperature. Such temperature 
independence at the early-time stage is a distinctive fea- 
ture that is not found in the simple Langevin dynamics. 
This indicates that dissipation at the early-time stage 
is dominated by a mechanism that is different from the 
temperature-dependent activation dynamics operating at 
the late-time stage. 

In order to help understand the temperature- 
dependence of the dynamic exponent, we have also exam- 
ined how ordering proceeds with time. Shown in Fig. |] 
are the snapshots at various time instants of ordering con- 
figurations in terms of the staggered chirality. Figure |] 
(a) and (b) exhibit the ordering process for quenching to 
low temperatures (T = 0.1) while (c) and (d) to higher 
temperatures (T = 0.4), where black and white regions 
represent domains with opposite (staggered) chirality or- 
dering. One can observe distinctive features in the do- 
main growth morphology in the two cases: At low tem- 
peratures, typical separation between consecutive corner 
defects along a domain wall grows larger in time. At 
higher temperatures, in contrast, the typical separation 
between corners remains more or less constant through- 
out the ordering process. Slow ordering at low tempera- 
tures (T = 0.1) compared with that at high temperatures 
(T = 0.4) can also be observed here, in agreement with 
the apparent temperature-dependence of l/z q obtained 



from the scaling analysis of the chirality correlation func- 
tion. 

The above morphological characteristics are related to 
the time dependence of the number N p of point defects 
and of the total length Ni of line defects. Simulations 
reveal power-law relaxation toward equilibrium with 
temperature-dependent exponents, as shown Fig. ^|(a) to 
(c). Namely, we have N p (t) = N p (t) — N p (oo) ~ t^"" and 
Ni(t) = Ni(t) - N t (oo) ~ t~ v \ where N p (oo) and A z (oo) 
are the values in equilibrium. These equilibrium values, 
which depend on the temperature, can easily be obtained 
from dynamic simulations, beginning with one of the the 
known ground state configurations. The resulting equi- 
librium values are negligibly small at low temperatures 
(T < 0.2) and the total amounts N p and Ni as well as 
the excess amounts N p and Ni display power-law decay. 
The corresponding decay exponents together with l/z q 
are shown in Fig. |5|(d) as functions of the temperature. 
At low temperatures, we have v p > v\ and point defects 
decay away faster than line defects, resulting in growing 
separation between corner defects along domain walls. At 
higher temperatures, on the other hand, we have v p « z/;, 
and the average distance between neighboring point de- 
fects along a line defect does not change -appreciably 
with time, similarly to the previous workstiULHI, which ar- 
gued, based on such numerical results, that there exists 
a roughening transition at a finite temperature. How- 
ever, it was argued in a recent work that due to the finite 
excitation energy for the double-step kinks, tho-domain 
walls should be rough at any finite temperature£3 If this 
is indeed the case, asymptotic decay exponents for the 
domain walls and corner charges should be the same at 
all finite temperatures. Unfortunately, numerical confir- 
mation of this would be very formidable due to the slow 
domain growth at low temperatures. 

V. SUMMARY 

We have studied the non-equilibrium relaxation and 
coarsening dynamics of two-dimensional superconduct- 
ing arrays in the overdamped limit under zero external 
driving current. In the case of unfrustrated arrays, we 
find that the friction constant of a vortex remains finite 
in the limit of large extent of the vortex; this is argued to 
lead to the absence of logarithmic factors in the length 
scale and in the vortex number decay, which is in rea- 
sonable agreement with simulation results. In the case of 
fully frustrated arrays, relaxation and coarsening dynam- 
ics are characterized by decay of line defects (chirality do- 
main walls) as well as of point defects (corner charges) . 
Here strong temperature dependence is found in the char- 
acteristics of the late-time domain growth, which can be 
explained by activated domain wall dynamics across the 
energy barriers with logarithmic size dependence. In the 
early-time regime, on the other handr, there exists a rel- 
atively long period of transient dynamics which is inde- 
pendent of the temperature; this may be attributed to 
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FIG. 5: Relaxation of the excess length of the line defects and of the excess number of point defects associated with staggered 
chirality domains in the array of size L — 128 at temperature (a) T = 0.20; (b) T = 0.30; (c) T — 0.40. Shown in (d) is the 
temperature-dependence of the decay exponents vi and v v for line defects and for point defects, respectively, together with the 
chirality growth exponent l/z q . 



the effects of initial random distributions of the strong dation through Grant No. 2000-015-DP0138. 
Lorentz force, making the thermal noise irrelevant. It 
would be of interest to investigate the coarsening dynam- 
ics of Josephson-j unction arrays under external magnetic 
fields and to compare with our simulation results. 
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